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C/3 ^ Abstract 

j^ , We present a simple and general argument showing that a class of dy- 

namical correlations give rise to the so-called Tsallis nonextensive statis- 
tics. An example of a system having such a dynamics is given, exhibiting 
i-rt , a- non-Boltzmann energy distribution. A relation with prethermalization 

r^ . processes is discussed. 

o 
o . 

1 Generalized statistics and correlations 

^ . In the recent paper by J. Berges et al. ^, the authors show that the dynamics 

CSJ ' of a quantum scalar field exhibits a prethermalization behavior, where the ther- 

C'i , niodynamical relations become valid long before the real thermal equilibrium 

is attained. There, although the system is spatially homogeneous, the single 
particle spectra are different from those determined by the usual Boltzmann 
^+ ' equilibrium. This means that the validity of thermodynamics does not imply 

f^ , thermal equilibrium in the sense of Boltzmann. 

Such a situation has also long been claimed by Tsallis |2l Bj , who intro- 
duced the g-entropy, defined in terms of the occupation probabilities {pa} of 
the microstate a by 

^: s,[{p^}]^ ^~^f\ (1) 



instead of the usual entropy S'[{pq}] = ^'^aPa^^Pa- AH thermodynamical 
relations are derived from the maximization of this g-cntropy. 



X\ 6S,-^5{E)^ + XSY,Po.^0, (2) 

a 

where 



{E}^ = %^ (3) 
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Figure 1: Schematic illustration of the effect of a strong correlation among two 
variables for minimizing a function. 



is the q-biased average of the energy. It is shown that, although the thermo- 
dynamical Legendre structure is maintained, the energy spectrum for such a 
system is not given by Boltzmann statistics but behaves asymptotically like a 
power in £^ |2inj. In other words, the thermal equilibrium in the sense of Boltz- 
mann is a sufficient condition, but not a necessary one for the thermodynamic 
relations to be valid. Extensive studies have been developed |3j and the ap- 
plications of Tsallis' entropy extend to biological systems, high-energy physics, 
and cosmology. However, in spite of these studies, the dynamical origin of the 
5-entropy and g-biased average of observables is not clear yet. 

It has been suggested that Tsallis' nonextensive statistics applies to a kind 
of metastable state which is realized before the true thermal equilibrium, or to 
some stationary state having long-range correlation |3|. This is exactly another 
important result of J. Berges et al., that their preequilibrium state is attained 
long before the real equilibrium is reached, that is, the system seems to at- 
tain first this preequilibrium, then slowly goes to the real equilibrium. This 
reminds us somewhat of the following situation which we encounter frequently 
in optimization problems. For example, let us consider the problem of finding 
a minimum of a function of two variables, as shown in Fig. 1. 

When the function has a narrow long valley, the usual steepest descent 
method always first goes to the closest point of the valley instead of going 
directly to the real minimum, as indicated by the dashed curve. The narrow 
valley means that there exists a strong correlation between the two variables. 
In other words, if there exist any strong correlations among variables, the mini- 
mizing path always tries to satisfy the correlation first, then follows downstream 
to the true minimum along the valley generated by the correlation. If the cor- 
relation is strong enough, non Boltzmann distribution may result as the stable 
configuration of the system (see Ref. .5J. 

The above consideration suggests that a preequilibrium state may correspond 
to a kind of narrow valley of the free energy of some subsystem generated by 
dynamical correlations of the system. If such a local minimum stays for a very 



long time, we may consider the system as in equilibrium, but the microscopic 
occupation probability is not that of the Maxwell-Boltzmann distribution. Then 
what characterizes such a local equilibrium? We suggest here that a general form 
of Tsallis' statistics is one answer. Of course, prethermalization involves a broad 
range of phenomena, and has been studied for a long time in several works (for 
recent references, see 0]). Here we concentrate on the discussion of how the 
states described by the Tsallis nonextensive statistics can be regarded as a kind 
of prethermalized states. 

Let us consider a system composed of N particles and let {Ea} denote the 
single-particle spectrum of the system and {Na} be the corresponding occupa- 
tion numbers. We further assume that the system is dilute enough that the 
interaction energies are negligible. Then, by definition, E = J^a ^aEa and 
A'' = J2a ^a ^^^^ respectively, the total energy and total number of particles of 
the system. Suppose there exists a strong correlation among any q particles of 
Na in any single-particle state a. Then the number of possible ways of forming 
such correlated subsystems in each a state is given by 

for Na ^ q. Note that here, q does not have to be integer. We call these 
correlated systems g-clusters. The total number of ways of forming g-cluster 
states for the whole system is 



M^J2^-'- (5) 



/ 77i\ Z^Q J'-ia^Q! L-iafa'^OL ,„•, 



If the g-clusters are formed with equal a priori probability, thus proportional 
to the number of ways of forming them, then the average energy per cluster 
particle is 

where Pa = Na/N is the probability of a particle occupying state a. 

Now, we further assume that the dynamical evolution of the system is such 
that, for a given number of g-clusters, the configuration of the system tries to 
minimize the energy of the correlated subsystem. That is, for some arbitrary 
initial condition, the system first generates correlations among particles in such 
a way to minimize the energy in the clusters. The narrow valley in Fig. 1 
corresponds to these "ground states" of the subsystem formed by g-clusters. 
Neglecting the component perpendicular to the narrow valley, the equilibrium 
corresponds to minima of the expectation value (E) , 

HE), = 0, (7) 



under the conditions 



a. 

M 



Y.PI 



m 



A = const. 



(8) 
(9) 



This is equivalent to the variational problem, 
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which determines the probabilities Pa- We obtain: 






l + ^(^a-(i?),) 



1/(1-9) 



(10) 



(11) 



where Zq = A^/'^'^^'^^ is a normalization factor. We see that this expression is 
nothing but the Tsallis distribution function, if we put A~^ = {q — I)/?- Note 
that there was no need to introduce the concept of entropy yet. The Tsallis 
entropy here corresponds to the simple normalization condition for the number 
of g-clusters. 

It is important to note that the minimization condition, Ea. HlU|) . of the 
energy of the correlated subsystems under the constraint Eq.(|5J) determines 
completely the occupation probabilities for the whole system as functions of A 
and of the single-particle spectrum, pa — Pa {A,{Ei}). When we change the 
external parameters of the system, such as the volume V , then the single-particle 
spectrum will suffer a change. However, if the dynamics of the system is such 
that the number of correlated particles q is kept constant during slow changes of 
V and N , then we can consider V , N, and A independent "state" parameters, 
since the occupation probabilities for the system are completely determined by 
these variables. We can thus calculate the total energy as function of them as 



E{V,A,N)^Nj2Po.Ea 



(12) 



Comparing the variation of this energy function with the First Law of Thermo- 
dynamics, 6E = —P6V + TSS + fiSN, we can identify the pressure 



— i:(i:(|£)/...^)^. 



the chemical potential /i = ^^PaEa, and the entropy term 



TdS 



"Edr)-.-" 



6A. 



(13) 



(14) 



Eq. H14fl . implies that we can relate the parameter A to the entropy of the 
system, up to an arbitrary function f{A), 



S - fiA) , 



and have a temperature given by 



T = 



fiA) 



E 



dA 



If we take J {A) = In A/ (1 - q), then 



Ea 



In(EaP^) 



(15) 



(16) 



(17) 



is the Renyi entropy. If instead we take f{A) = (A — 1)/(1 — q), S corre- 
sponds exactly to the Tsallis entropy (^3). This indeterminacy is obvious also 
from Eq. UlUfl . since the constraint term A(^p^ ~ -A) can also be taken as 
•^' [/ (SPq) ~ /(^)] without changing its significance, for any single-valued 
function /. Thus, it is clear that the exact form of / is not essential in our 
approach. 

We should stress that the thermodynamical properties of the preequilibrium 
state discussed above are not the nonextensive ones discussed by Tsallis et al. 
12 131 . There, the thermodynamical and statistical properties for g-modified 
values are studied. In the present approach, the g-modified quantities refer 
to the subsystem of g-clusters. When we want to calculate, for example, the 
internal energy of the system, the relevant quantity is not {E) but just the usual 
total energy E. Therefore, the effective equation of state, that is, the relation 
between pressure and energy density given by Eq. (|13|l . will apply immediately 
as the real physical relation which can be used, for example, in hydrodynamical 
calculations. However, the "entropy" here is not the real entropy of the whole 
system but just refiects the total number of clusters among correlated particles. 
Thus, the time evolution of this quantity is not guaranteed to obey the Second 
Law of Thermodynamics, that is, dS/dt > 0. Hereafter, in order to distinguish 
from the real entropy, we call this S the "pseudo-entropy" of the total system. 

On the other hand, note that the variational equation Eq. I|l()(l can also be 
obtained as that of maximizing the cluster multiplicity M at a fixed average 
energy per cluster. We would have the variational equation. 



M 



■(3{{E), 



Er 



J2p° 



0. 



(18) 



where Eq is the average energy per cluster, and (3 and v are Lagrange multipli- 
ers. This leads to the same equation as Eq. (|10|l . Since M is the total number 
of possible ways to form clusters, we can define the "entropy" of the system 
of clusters by Sq = In M, and this corresponds to the Renyi entropy, Eq. ((T7|) . 
Therefore, Eq. Hll)|) is nothing but the maximization of the entropy defined only 



for the restricted subspace corresponding to the degrees of freedom of clusters 
in the entire phase space of the system. 

In the above discussion, we considered the value of q, i.e., the number of 
strongly correlated particles as constant for any single-particle state a. However, 
in a more general case, they may be different for each single-particle state and 
may even depend on the proper occupation numbers Pa- Therefore, for a more 
general form of the correlation pattern, the number of cluster states M^ should 
be written as 

Af„ = M„ [{p^ }] = /" (19) 

where Fa [{Pa}] are functionals of {pa} ■ Then, the Eq. H1U|) should have the 
form 



EMa[{Pc.}]Ea 



EMaiiPc}] 



+ \(J2 ^^" [{p"}] - ^^) +" {J2p-' " ^)] = • (20) 



The resulting probabilities Pa will not be of the form of Eq. (|11() anymore, yet 
functions of V, N, and M. This suggests that the preequilibrium state can 
be attained for a wide variety of single-particle spectra and the corresponding 
pseudo-entropy will be given as a function of M, S — f (M) — f {J2a^oi)- 
This is presumably the most general form of the pseudo-entropy. However, 
we should recall that it does not necessarily possess the properties of entropy 
required by the usual thermodynamics because the state specified by Ea.lj^ 
does not correspond to the conventional equilibrium of the system. 

2 Example 

Discussion in the previous section indicates what kind of dynamical correlations 
may lead to Tsallis type occupation probability distribution. Just to see a 
concrete example let us consider a very simple toy model described below. 

1. Initially TV particles are distributed as {No,Ni, ...,Na, ...} over equally 
spaced energy levels. Let us denote the energy of the i— th particle as Ei. 

2. Choose randomly a pair of particles, say, i and j. 

3. Energies of i and j are updated according to one of the following alterna- 
tives: 

(a) The new energies are set to the lower one of Ei and Ej, that is, 
Ei,Ej ^> E[ — Ej — Min{Ei,Ej), then choose another k-th{k ^ i,j) 
particle randomly and attribute E'/. = Ek+Max {Ei,Ej)—Min{Ei, Ej) 
to conserve the total energy. 

(b) Change the energies as E^ — > E'^ = E^ ± AE and Ej -^ Ej = Ej =F 
AE, where AE is the level spacing. Here, if one of El,E', becomes 
negative, then this step is skipped. That is, we have to keep always 
E, > 0. 



4. The alternatives (a) and (6) are chosen randomly, but the ratio of the 
average frequency of (a) to (b) is kept as a constant, r. 

It is well-known that for r = in the above model, the ultimate single- 
particle spectrum will be the Boltzmann distribution. When r > 0, the collisions 
of type (a) are forming a kind of cluster, while the collisions of type (6) may 
destroy a correlation formed before. Thus, we expect the number of correlated 
particles will reach some stationary value. Furthermore, by construction, the 
energy of a correlated pair always tends to diminish. In this way, we may expect 
that the above system will lead to the situation described by Eq. 11U|) . In Fig. 2, 
we show the results of simulations, for several values of r from 0.00001 to 0.1. It 
is interesting to note that r > leads to a non-Boltzmann distribution which is 
well approximated by the Tsallis distribution, pa = C/ [1 + {q — 1) (3Ea] ^ , 
where /3 = 1/(3 - 2q) (E) and C = (2 - q)/(3 - 2g) (E) are determined from the 
normalization condition and conservation of energy. One parameter fits with 
respect to q were performed and the results are indicated by the continuous 
curves in this figure. For r — *■ 0, the spectrum converges to the Boltzmann dis- 
tribution Pa -^ e^^°/^^^ as expected. Note that for a one dimensional case like 
the present model, the Tsallis distribution is valid only for q < 3/2, otherwise 
the energy expectation value diverges. For larger values of r, the fitted value of q 
tends to this limiting value, but the distribution begins to deviate substantially 
from the Tsallis distribution in the low energy region. 

It is important to note that these distributions are the stationary and stable 
ones. We confirm that starting from any different initial conditions, the system 
always converges to the same final distribution, uniquely determined by the r 
parameter. Furthermore, the convergence for r > is much faster than for 
r = 0. For example, taking r = 0.01 we find that the distribution converges 10 
times faster than for r = 0. 

In the above toy model, time reversal is violated in type (a) collisions. How- 
ever, this is not the crucial factor to obtain the non-Boltzmann distribution. We 
have checked this in a more elaborate model which has time-reversal invariance. 

3 Conclusions and perspectives 

We have shown that a very simple mechanism of dynamical correlations among 
particles may lead to a state described by the so-called nonextensive statistics 
of Tsallis. It is shown that the nonextensive statistics can be interpreted as the 
statistical mechanics for a subspace defined by these strongly correlated clusters 
in the entire phase space of the whole system. Furthermore, the requirement 
of the minimum energy per cluster or the maximum "entropy" of the subspace 
defines the thermodynamical properties of the whole system too. There, the 
energy is defined by the usual sum of single-particle energies multiplied by their 
occupation probability, and not by the g-weighted average of the energy. We 
have also shown a toy model for which such a mechanism actually takes place. 
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Figure 2: Energy spectrum after a large number of collisions per particle, start- 
ing from a distribution peaked at (E) — 22.5. The probability of cluster-forming 
collisions ranges from r = 0.1 (top) to r = 0.00001 (bottom). Tsallis distribu- 
tions fitted to the calculated spectra are also shown, together with the corre- 
sponding g-values. 



In this example, the stationary state is attained with a single-particle spec- 
trum which is well represented by that of Tsallis statistics. These stationary 
non-Boltzmann distribution are attained much earlier than the case of no cor- 
relations, hence the Boltzmann equilibrium. The reason for reaching the sta- 
tionary distribution much earlier than the real thermal equilibrium is that the 
formation of correlated clusters serves as a kind of catalyzer to distribute the 
energy among particles more efficiently, involving more than 2-body collisions. 
In this aspect, the dynamical mechanism responsible for the non-Boltzmann 
distribution is similar to that of the kinematical approach discussed by Rafel- 
ski et al. [H] . Many examples of systems for which Tsallis statistics have been 
applied may also be interpreted in our scheme easily. We also found that this 
non-Boltzmann equilibrium may happen not only with a power type but with 
a much wider class of single-particle spectra. 

The present work is rather phenomenological and so the conclusion is still 
speculative. More theoretical work based on realistic field theoretical models 
should be developed. Further investigations in this line are in progress. This 
work has been supported by CNPq, FINEP, FAPERJ, and CAPES/PROBRAL. 
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